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Abstract. The latest version of the planetary ephemerides developed at the Paris Observatory and at the Besangon 
Observatory is presented here. INPOP08 is a 4-dimension ephemeris since it provides to users positions and 
velocities of planets and the relation between TT and TDB. Investigations leading to improve the modeling of 
asteroids are described as well as the new sets of observations used for the fit of INPOP08. New observations 
provided by the European Space Agency (ESA) deduced from the tracking of the Mars Express (MEX) and Venus 
Express (VEX) missions are presented as well as the normal point deduced from the Cassini mission. We show 
the huge impact brought by these observations in the fit of INPOP08, especially in terms of Venus, Saturn and 
Earth-Moon barycenter orbits. 

Key words, celestial mechanics - ephemerides 



1. Introduction 

Since the first release, INPOP06, of the plane- 
tary ephemerides developed at Paris and Besangon 
Observatories, (Fienga et al. 2008, www.imcce.fr/inpop), 
several improvements have been conducted on the dynam- 
ical modeling of the INPOP ephemeris. The observation 
dataset has also been substantially increased, especially 
with the addition of ranging data from the ESA space 
missions Mars Express and Venus Express. The resulting 
new version of INPOP planetary ephemerides, INPOP08, 
is presented here with the description of its new features. 
The INPOP08 ephemeris has also been fitted to all avail- 
able Lunar Laser Ranging data (Manche et al., 2007). 

One of the novelties present in INPOP08 is the ad- 
dition, in the distribution of the ephemeris, of a time 
scale transformation TT-TDB that is coherent with the 
ephemeris. The basic idea is to provide to users po- 
sitions and velocities of Solar System celestial objects, 
and also, time ephemerides relating the Terrestrial time- 
scale, TT, and the time argument of INPOP, the so-called 
Barycentric Dynamical Time TDB, based on the defini- 
tion adopted by the International Astronomical Union in 
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2006. Such a release of planet and time ephemerides enable 
us to go towards four-dimensional planetary ephemerides. 
Section [2] describes the INPOP procedure for the compu- 
tation of the TT-TDB relation. 

Section [3] is devoted to a brief account of the new con- 
straints and modeling implemented in INPOP for the as- 
teroid perturbations. As in INPOP06, 300 asteroids are 
included in the dynamical equations of INPOP08, and the 
remaining ones are modelized as a ring. With respect to 
INPOP06, the ring model and the asteroid selection have 
been improved, and the precise description of these ad- 
vances are made in (Kuchynka et al, 2009). 

An important part of the new observations used for 
the INPOP08 fit consists in the tracking data provided by 
the ESA space missions Mars Express and Venus Express 
(Morley 2006a, 2007a, 2007b). These datasets are the first 
radio ranging data provided by ESA, and their acquisition 
and reduction process is described in section [H VEX ob- 
servations bring a very important set of informations on 
Venus orbit. It is especially of interest since this orbit is 
much less perturbed by asteroids than Mars and therefore 
has greater potential for precise ephemeris motion, funda- 
mental physics testing and reference frame establishment. 
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The section 5 deals with the INPOP fit obtained by 
comparisons between the dynamical modeling and the ob- 
servations. In this process, 34 asteroid masses were fit- 
ted against only 5 in INPOP06. Comparisons are made 
between obtained asteroid masses and other published 
masses. Values for the fitted Earth-Moon barycenter and 
Sun oblateness J2 are also given. We have also fitted the 
AU and comparisons are provided with the latest deter- 
minations of DE414 (Standish 2006; Konopliv et al. 2006) 
and DE421 (Folkner 2008). 

In addition, we have performed as well a second fit, 
where the AU is given the lERS Conventions 2003 value 
(IERS03), and the Solar CM is deduced. 

On several occasions, for planned high precision ob- 
servations, we had some enquiries about the real accu- 
racy of the position or velocities given by the plane- 
tary ephemerides. This is actually a difficult question 
to answer, but we have tried in section [6] to provide 
some estimates of these uncertainties, by comparison with 
INPOP06 and DE421. The last section is devoted to the 
conclusions and perspectives. 



2. INPOP as a 4-D planetary ephemeris 

With INPOP08, we aim to produce planetary ephemerides 
as fully compatible as possible with the relativistic back- 
ground recently adopted by the astronomical community 
and summarized by the IAU2000 and IAU2006 conven- 
tions (Soffel el al., 2003). This leads to the production of 
an ephemeris in TDB, and to the construction of a TT- 
TDB transformation. The following subsections describe 
the various steps involved in this process and the impact 
for the ephemeris users. 



2.1. The TCB-TCG transformation 

Two reference systems are defined: a global one, the 
Barycentric Celestial Reference System (BCRS), covering 
the whole Solar System and a local one, the Geocentric 
Celestial Reference System (GCRS), which is physically 
suitable for the modeling of processes in the vicinity of 
the Earth. BCRS is particularly useful when one wants to 
model the light propagation or motion of celestial objects 
in the Solar System. We can then, in the mass-monopoles 
approximation, write the equation of motion of bodies 
as well as the conservation laws satisfied by the Solar 
System barycenter (see Damour and Vokrouhlicky 1995); 
all these features being already implemented in INPOP06. 
(Damour and Vokrouhlicky 1995) gives the complete con- 
servation equations based on (Damour, Soffel and Xu 1991) 
formalism including spin/spin, mass/mass and mass/spin 
couplings. However, each fundamental reference system 
has its own time-scale: TCB for the BCRS and TCG for 
the GCRS. The relation between TCB and TCG can be 
derived, but time transformations are only determined for 
specified space-time events (one coordinate time and three 



spatial positions) so, at the geocenter, the transformation 
reads (Damour, Soffel and Xu 1991, Soffel et al. 2003): 

dTCC 1 1 / 1 \ 

^-^^-.<^iTCB)^-,PiTCB)^o(^-^) (1) 

c being the speed of light in a vacuum and where 
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Capital latin subscripts A and B enumerate massive bod- 
ies, E corresponds to the Earth, AIa is the mass of body 
A rEA = xe - xa, Tea = \rEA\, va = xa, aA = va, 
a dot signifying time derivative with respect to TCB and 
XA being the BCRS position of body A. 

Eqs. ([!]), ^ and ^ exhibit that the time transfor- 
mation between TCB and TCG is explicitly related to 
the positions and velocities of Solar System bodies, so 
to the planetary ephemeris itself. As stressed by Klioner 
(2008), if we are distributing a time transformation to- 
gether with positions and velocities, we are building a 
four-dimensional planetary ephemeris. This time transfor- 
mation will be called in the following time ephemeris to be 
close to the terminology of Irwin and Fukushima (1999). 

An issue remains. The use of TCB in planetary 
ephemerides should induce important changes in numer- 
ical values of planet masses, initial conditions and astro- 
nomical unit commonly adopted by users. 

Therefore, even if a TCB-based ephemeris and a TDB- 
based ephemeris are related linearly, caution has to be 
taken when one wants to shift from one ephemeris to an 
other as it will induce large changes for all users. For a 
more complete discussion, see for instance (Klioner 2008). 

2.2. The TT-TCG transformation 

The lAU time realization is done by the International 
Atomic Time (TAI) with a service running since 1958 and 
attempting to match the rate of proper time on the geoid 
by using an ensemble of atomic clocks spread over the 
surface and low orbital space of the Earth. TAI is con- 
ventionally related to the Terrestrial Time (TT) and the 
Universal Coordinate Time (UTC) by the two following 
relations: 



TT{TAI) = TAI + 32.184s . 



(4) 
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and 



2.4. Implementation in INPOP08 



UTC = TAI - leap seconds , 



(5) From equations ([T]), (O and (O, we have (Klioner, 2007): 



leap seconds being added by the International Earth 
Rotation Service at irregular intervals to compensate for 
the Earth's rotation irregularities. Mainly, leap seconds 
are used to allow UTC to closely track UTl, which effec- 
tively represents the Earth rotation. 

It is thus more convenient to consider a regular time-scale 
like TT. Moreover due to its IAU1991 definition, TT is 
related to TCG by a fixed linear function as follows: 



TCG -TT^ LgX {JD - 2443144.5) x 86400 , 



(6) 



where Lg = 6.969290134 x 10"^° and JD is TAI measured 
in Julian days. 

2.3. The TDB-TCB transformation 

The TDB has a more tumultuous history. It was intro- 
duced by IAU1976 in order to remain close to TT up to 
periodic variations. However such a definition was flawed 
because in that case TDB cannot be a linear function 
of TCB. Consequently, the relativistic equations derived 
in TCB cannot be simply adapted to TDB. In practice, 
Caltech/JPL, Harvard/CfA, and RAS/IAA ephemeris 
programs have used a relativistic coordinate time whose 
mean rate is automatically adjusted to the mean rate of 
TT by way of the ephemeris fits. As such, those times, 
often referred to as Tep/i, can be related to TCB by the 
following equation with the notation of equation [TJ 

dTeph . 1 



TCB 



= 1 + -a{TCB) 



and then from (Irwin and Fukushima 1999), 



TCB = 



1 



(7) 



(8) 



Differences between T^ph and TCB have been estimated 
by numerical quadrature method of equation [7] which 
is equivalent to equation [TUl limited at the terms. 
Chebychev polynomials of the solutions were provided to 
users and corrections were brought to the analytical repre- 
sentation of TDB done few years earlier by Fairhead and 
Bretagnon (1990). 

This situation changed with the recommendation B3 
of IAU2006 which made TDB a fixed linear function of 
TCB as follows: 



TDB = TCB - Lb{JD - To) x 86400 + TDBq, 



(9) 



where To = 2443144.5003725, Lb = 1.550519768 x 10-^ 
TDBo = -6.55 x lO'^s and JD is the TCB Juhan date. 
TCB value (like TT and TCG ones) is Tq for the event 
1977 January 1 OOh 00m 00s TAI at the geocenter. It in- 
creases by one for each 86400s of TCB. The use of a TDB 
planetary ephemeris as defined in equation ([9]) has no sig- 
nificant impact for normal users: values of masses and ini- 
tial conditions of solar system objects being the same at 
the common level of accuracy. 



d {TT - TDB) 
dTDB 



= [LB + -^a] {1 + Lb-Lg)~ Lg 



(10) 



One can notice that the values of a and /? do not 
change when using quantities (CM, positions, velocities 
and accelerations of bodies) expressed in TDB instead 
of TCB units. It is then straightforward to construct a 
time ephemeris TT - TDB = /(TDB) (noted ATDB = 
f{TDB)). This is the choice done for INPOP08. 

The central part of our implementation into INPOP is 
to numerically integrate Eq. pO]) together with the equa- 
tions of motions of all bodies. Because the right member of 
(fTO]) does not depend on ATDB, the equation (fTO|) is not 
strictly an ordinary differential equation. In the compu- 
tation of a, A enumerates all massive bodies of the Solar 
System, that is, the Sun, the planets, Pluto, the Moon and 
all the 303 asteroids. In /3 (which is divided by c"*, and so is 
less important than a), A and B enumerate all bodies ex- 
cept the 298 "small" asteroids. For the acceleration term 
a A, all the interactions are taken into account, including 
newtonian interactions, relativistic corrections, figure and 
tide effects. It was already needed for the equation of mo- 
tion of the corresponding body, and no additional work is 
therefore necessary to compute it. 

Because the (TT - TDB) value is unknown at J2000 
(it depends on the ephemeris), the initial condition is set 
to zero. The quantity integrated in the state vector (in- 
cluding positions and velocities vectors of bodies) is then 
ATDB + k where k is an offset determined later, just be- 
fore building the Chebychev polynomials, by using the fol- 
lowing condition: for the event 1977 January 1st OOh 00m 
00s TAI at the geocenter, TDB julian day is To + TDBq 
and ATDB = -TDBq. 

At this point, the difference TT — TDB can be com- 
puted from any value of TDB. But in the reduction pro- 
cess of observations and because they are dated in UTC 
time scale, the transformation TT — TDB as a function 
of TT (noted ATT) is needed. 

A similar differential equation as eq. (jlOp can be found 
for d{ATT)/dTT (Klioner, 2007), but it is not neces- 
sary to integrate it. From the relation (TT — TDB) = 
f{TDB) computed with INPOP, one can notice that 
(TT - TDB) = f{TT - (TT - TDB)). ATT is then solu- 
tion of an implicit equation, which can be solved by itera- 
tions. In fact, only one is necessary, discrepancies between 
TT and TDB being smaller than 2 milliseconds. 

The formal differences between the present procedure 
and TE405 (Irwin & Fukushima, 1999) are small as il- 
lustrated in Fig. [T] The linear drift is measured to 6.74 
nanosecond per century (ns/cy). This value is consistent 
with the IAU2006 resolution B3 (less than 1 nanosecond 
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Fig. 1. Differences (in nanosecond) between TE405 (cor- 
rected from the 6.55 x 10"^ second offset) and the (TT- 
TDB) integrated in fNPOPOS. X axis is Terrestrial Time, 
expressed in years from J2000. 

per year between TT and TDB). These small discrepan- 
cies could be surprising because TE405 does not take into 
account the terms in l/c"*. Neglecting them in eq. PTI)1 
induces an important drift of 346.0 ns/cy , close to the 
value AL^^-* = 346.2 ns/cy from (Irwin & Fukushima, 
1999). But this drift can be compensated by a change of 
the constant is or Lq- Now that TDB is defined as a con- 
ventional and fixed linear function of TCB (see eq. O and 
Lb is fixed to the value given by the lAU 2006 resolution 
B3, terms in 1/c'* are thus essential. 

Because discrepancies are small, no impact is no- 
ticeable in the differences between observed planet posi- 
tions and positions computed with TE405, the relation of 
Fairhead and Bretagnon (1990), or with the Chebychev 
polynomials representing the INPOP TT^TDB. No ad- 
ditional iteration is then necessary, even if at each step of 
the adjustment of INPOP to observations, the computa- 
tion of TT^TDB coefficients is automatically iterated. 

3. New constraint on asteroid modeling 

3.1. Supplementary selection of asteroids 

For INPOPOS, we also slightly revised the INPOP06 se- 
lection of asteroids perturbing Mars and the inner plan- 
ets. This was done with the same method as the estima- 
tion of the mass of the ring described in the [2] section : 
24635 selected asteroids are assigned with reasonable dis- 
tribution of masses according to available data and the 
Statistical Asteroid Model (Tedesco et al. 2003). A Monte 
Carlo study allows to assign to each asteroid the probabil- 
ity of being among the 300 most perturbing asteroids in 
terms of amplitude of the perturbation on the Earth-Mars 
distance between 1969 and 2010: For each of the 24635 as- 
teroids, two integrations have been done between 1969 and 
2010: one with the asteroid i and one without the asteroid 
i, i varying from 1 to 24635. The differences between the 
two integrations give the impact induced by the asteroid 
i. To built our list of asteroids, we have studied the impact 
not only on the Earth-Mars distances, but also the Earth- 
Venus and Earth-Mercury distances. We are thus able to 
compile the most probable list of the 300 most perturbing 
asteroids and compare it to the INPOP06 selection used 



m INPOPOdll. Three asteroids have predicted perturba- 
tion amplitudes well over 30m and are absent from the 
INPOP06 selection: 60 Echo (amplitude reaches 170m), 
585 Bilkis and 516 Amherstia. We added these to the 300 
already integrated in INPOP06 and thus there are in total 
303 asteroids integrated individually in INPOPOS. More 
details on the methodology used for obtaining the most 
probable list of the top 300 perturbing asteroids can be 
found in (Kuchynka et al. 2009). 

3.2. The ring 

In the former version of INPOP (Fienga et al. 2008), per- 
turbations of asteroids on planets were modeled with 300 
individual asteroids and a static circular ring at 2.8 AU. 
Five asteroid masses, 3 taxonomic densities (attributed to 
the remaining 295 asteroids) and the mass of the ring were 
fitted to observations. 

After the integration of such a model on a 100-years 
time interval, it appeared that the asteroid ring induced a 
drift of several meters in the position of the Solar System 
barycenter (Figl2]). The static ring was then replaced by 
a more realistic implementation that now conserves the 
total linear and angular momenta of the system: the ring 
interacts fully with the planets and is no longer assumed 
to act only in their ecliptic planes. Its center is attached 
to the Sun and its orientation is an integrated parameter 
which evolves with time. The ring's interaction with other 
objects is actually identical to the interactions of an aster- 
oid on a circular orbit averaged over the asteroid's mean 
orbital motion. Besides eliminating the barycenter drift, 
the advantage of the new implementation is that the ring 
is taken into account in a more realistic way and its pres- 
ence in the model is more meaningful. The ring's radius 
was chosen at 3.14 AU. 

Another novelty is that the mass of the ring is not 
fitted to observations but estimated independently. This 
estimation is made by calculating the amplitude of the 
perturbation on the Earth-Mars distance exerted by all 
main-belt asteroids but the 303 most perturbing ones de- 
termined previously. 24635 asteroids are considered as a 
model of the main belt and a simple scheme based on 
the Statistical Asteroid Model (Tedesco et al. 2005) is 
used to assign each asteroid with a reasonable distribu- 
tion of masses. A Monte Carlo experiment where aster- 
oids are assigned random (but reasonable) masses allows 
to calculate the corresponding perturbation from all the 
asteroids but the 303 most perturbing ones. For each 
random set of masses, a ring's mass is determined so 
as to fit the perturbation of the ring to the global ef- 
fect. This leads to the estimation of the ring's mass at 
M„„g = (1 ± 0.3) X W-^°Mq for a ring at 3.14 AU 
(Kuchynka et al. 2009). 



^ the list of asteroids integrated individually in INPOP06 or 
DE405 is given in (Standish and Williams 2001) 
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Time from J2000 (year) 

Fig. 2. Impact of the asteroid ring on the position of the 
SSB over 100-years integration. 



4. Presentation of MEX and VEX observations 

The tracking of Mars Express (MEX) and Venus Express 
(VEX), while orbiting their planetary namesakes, com- 
prises two-way coherent Doppler and range measurements. 
The Doppler data are converted to range-rate: the com- 
ponent of the spacecraft's velocity relative to the ground 
station along the direction from the station to the space- 
craft. The determination of the spacecraft orbits uses only 
the range-rate data since, under normal circumstances, 
the additional inclusion of the range data leads to only 
insignificant improvement in the accuracy of the orbit so- 
lutions. 

Within the orbit determination, the computed val- 
ues of the observables rely on high fidelity modeling of 
the dynamics and the signal path (Budnik et al. 2004; 
Moyer 2000). As part of this modeling, the orbital states 
of Mars and Venus are taken from the JPL DE405 plan- 
etary ephemerides (Standish 1998b). Over time the range 
residuals exhibit a signature whose predominant contri- 
bution is the error in the distance from the Earth to the 
planet as given by the planetary ephemerides. The range 
residuals are therefore derived data that are very useful 
for improving the accuracy of planetary ephemerides. The 
accuracy of the range residuals depends on a number of 
factors, which are described in the following subsections. 

4.1. Accuracy of the MEX orbit determinations 

MEX, the first European mission to Mars, was launched 
on 2nd June 2003 and inserted into Mars orbit on 25 
December 2003 and is presently extended until May 2009. 

MEX is tracked with both ESA 35 m deep space anten- 
nas and those of the NASA Deep Space Network (DSN) 
(mainly at the Goldstone complex), in almost equal pro- 
portions in terms of tracking duration. Most of the ESA 
tracking is from New Norcia (NNO), in Western Australia, 
the rest from Cebreros, near Madrid in Spain. 



The quantity of DSN range data and hence range resid- 
uals is far higher than the quantity of ESA range residuals 
due to the extraction of one ESA raw data point every 20 
minutes whereas DSN range measurements of MEX are 
made at intervals of 3 minutes 27 seconds. The standard 
deviation of range residuals for individual passes gives a 
good indication of the random noise on the measurements. 
For ESA data, the average standard deviation of the (two- 
way) residuals is less than 1 m. For NASA/DSN data it is 
below 0.5 m. 

In terms of orbit geometry, the operational orbit is 
near polar and elliptic with an apoapsis altitude of a lit- 
tle over 10000 km and a periapsis altitude that has varied 
between 250 km and 340 km. The orbital period averaged 
6.72 hours until late in 2007 when a series of five ma- 
noeuvres at periapsis increased the period to 6.84 hours. 
MEX is equipped with a fixed high-gain antenna (HGA) 
that must be pointed towards the Earth during tracking 
passes and tracking is done in X-band up- and downlink. 
The primary science data collection is around periapsis 
passage when the spacecraft HGA is not Earth-pointing. 
For orbit determination purposes, the information content 
within the Doppler data is highest around periapsis, so the 
orbit solution accuracy is adversely affected by the lack of 
such data. The orbit determination accuracy is also lim- 
ited by other factors; most notably the imperfect calibra- 
tion of perturbing velocity increments caused by thrust- 
ing, to off-load the accumulated angular momentum of the 
reaction wheels, and also the problems in accurately mod- 
eling small forces due to solar radiation pressure and the 
tiny but highly variable atmospheric drag at each periap- 
sis passage. 

For much of the MEX operational mission, routine or- 
bit determinations were made twice per week based on 
tracking data arcs of 5-7 days duration, corresponding 
to approximately 18-25 orbital revolutions, with a typi- 
cal overlap of 2 days between successive arcs. Nowadays, 
orbit determination is made weekly using an arc of about 
10 days duration and therefore with similar overlaps. The 
differences between the values of common range residuals 
computed from successive orbit determinations provide an 
indication of the effect of orbit determination errors on the 
accuracy of the residuals. Outside of periods of superior 
solar conjunction, these differences are almost always be- 
low 3 m. The range residuals that are retained are taken 
from the earlier of each overlapping solution. 

4.2. Accuracy of the VEX orbit determinations 

VEX, the first European mission to Venus, was launched 
on 9th November 2005 and inserted into Venus orbit on 
11th April 2006. Its mission is presently extended until 
May 2009. 

The primary ground station supporting VEX is 
Cebreros, which tracks the spacecraft almost every day. 
Between 27th April 2006 and 13th March 2008, there were 
663 Cebreros passes whose accumulated duration accounts 
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for 97.7 % of the total tracking duration. New Norcia 
has provided two-way coherent radiometric data during 
12 passes and NASA DSN stations during 26 passes, of 
which 16 were made with DSS 43 at Canberra, mainly 
during the superior solar conjunction in autumn 2006. The 
sampling rates of the range data and the random noise on 
the measurements are the same as for MEX. 

The operational orbit is polar and highly elliptic. 
There is hardly any precession of the line of apsides, so 
the periapsis, whose argument is currently at 95°, remains 
close to the planet north pole. The apoapsis altitude is 
about 66500 km and the periapsis altitude has been con- 
trolled to stay within the range of 185 km to 390 km. The 
orbital period is nominally 24 hours and the maximum 
excursion has never been more than 6 minutes. Every or- 
bital revolution, in the vicinity of apoapsis but not within 
ground station visibility, the momentum of the reaction 
wheels is off-loaded by thrusting. The perturbing AV into 
the orbit is typically in the range 15-25 mm/s which is 
substantially higher than for the MEX wheel off- loadings. 
The spacecraft attitude and the direction of the thrust 
are chosen so that the manoeuvres help to control the or- 
bit phasing. The control is such that the signal elevation 
from the daily Cebreros passes rises to 10 ° (when telecom- 
manding may start) always close to 2 hours after periap- 
sis passage. The pass ends either 10 hours later, close to 
apoapsis, or at 10° descending elevation, whichever is ear- 
lier. Tracking data obtained in X-band up and downlink 
are thus almost never acquired during the descending leg 
of the orbit, nor around periapsis. The combination of the 
unfavourable pattern of tracking data arcs, imperfect cal- 
ibration of the wheel off-loadings and deficiencies in mod- 
eling forces due to solar radiation pressure, together with 
other factors depending on the nature of the orbit cause 
the accuracy of the orbit determination to be worse than 
that for MEX. Typical values of the differences between 
range residuals derived from successive orbit solutions are 
a few metres but, occasionally, even away from solar con- 
junction periods, the differences can reach as high as 10 
m. 

4.3. Spacecraft transponder group delay 

Subtracted from each range measurement is the nominal 
value of the group delay of the on-board transponder. For 
MEX, the value corresponding to the normally used X- 
band up- and downlink signals is 2076 nanoseconds (about 
622 ni) . For VEX, that has a virtually identical transpon- 
der, the value is 2085 ns (about 625 m). The nominal 
value is the average value measured at different occasions 
on ground before launch. From the variations in measured 
values, ostensibly made under identical conditions, it is 
thought that the systematic error of the nominal value 
should not be larger than 30 ns (about 10 m). In addition, 
it is known that the group delay is not perfectly stable and 
can fluctuate by a few ns, depending upon variations in 
a number of parameters such as temperature and signal 



strength. For MEX, these error estimates appear reason- 
able, and perhaps a little conservative, based upon the 
consistency with results determined from the NASA MGS 
and MO spacecraft range data during the spring of 2005. 
For VEX, there is no independent means to verify the er- 
ror estimates. 

4.4. Superior solar conjunction 

During the time that MEX and VEX range residuals have 
been generated and archived, both missions have experi- 
enced a superior solar conjunction. The MEX Sun-Earth- 
probe (SEP) angle remained below 10° for two months 
centred on 23rd October 2006, when the minimum SEP 
angle was 0.39° (1.6 solar radii). The VEX SEP angle 
was continuously less than 8° over two months centred on 
27th October 2006, when the minimum SEP angle was 
0.95° (3.8 solar radii). The effects on spacecraft radio- 
metric data at these conjunctions have been described by 
(Morley and Budnik, 2007). The signals to and from the 
spacecraft pass through the solar corona surrounding the 
Sun. The free electrons in the plasma cause a group delay 
on ranging measurements. Since the electron density in- 
creases with decreasing distance from the Sun, following, 
at least approximately, an inverse square law, the delay 
increases as the SEP angle diminishes. No solar corona 
model was applied when computing the range residuals, 
so the increased delay is the cause of the peak in figure [3] 
and gaps in figured The existing solar corona models that 
could be used to correct the range residuals are not very 
accurate. They cannot take into account the quite large 
day-to-day variations in the signal delay caused by short- 
term fluctuations in solar activity like sunspot formation, 
flares and coronal mass ejections, all of which can influ- 
ence the surrounding electron density. A secondary cause 
of increased errors in range residuals at solar conjunction 
is due to the main effect on Doppler measurements of a 
substantial increase in noise. When the SEP angle falls to 
about 1° , the measurement noise typically increases up to 
two orders of magnitude higher than is usual at large SEP 
angles. The accuracy of range residuals is then indirectly 
and adversely affected by a degradation in the accuracy of 
the orbit determination solutions. That is why the deci- 
sion was taken to omit 30 days before and after the solar 
conjunction in range residuals. To illustrate the impact of 
solar conjunction, peak in residuals is plotted in the case 
of VEX data on figureO On figurelH the gap in MEX data 
around October 2006 is induced by the solar conjunction. 
Such solar conjunctions also occured during the MGS and 
MO missions and caused observation gaps such as about 
January 2000 and mid-2002. 

5. INPOP08 Fit 

For the fit of INPOP08, MEX and VEX observations pro- 
vided by ESA and described previously were added to the 
INPOP06 data sets used for its adjustment. See (Fienga 
et al. 2008) for a more detailed description of this data 
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Table 1. Residuals obtained from INPOP06 and 
INPOP08 



Planet 


Data 


Nbr 


INPOP06 


INPOP08 




[unit] 




Icr 


Icr 


Mercury 










Direct radar 


range [m] 


462 


894 


842 


Venus 










Magellan 


VLBI [mas] 


18 


2 


2 


Direct radar 


range [m] 


488 


1386 


1376 


VEX 


range [m] 


15131 


185.829 


4.6 


Mars 










MGS/MO 


range [m] 


10410 


5.954 


1.57 


MEX 


range [m] 


6006 


13.56 


2.07 


Path 


range [m] 


90 


7.63 


12.46 


Vkg 


range [m] 


1245 


17.396 


18.53 


Mixed 


VLBI [mas] 


96 


0.4 


0.4 


Jupiter 










Galileo 


VLBI [mas] 


24 


12 


11 


Optical 


ra [arcsec] 


5616 


0.343 


0.343 


Optical 


de [arcsec] 


5534 


0.332 


0.338 


Saturn 










Optical 


ra [arcsec] 


5598 


0.347 


0.346 


Optical 


de [arcsec] 


5573 


0.312 


0.311 


Cassini 


ra [mas] 


31 


5 


4 


Cassini 


de [mas] 


31 


7 


7 


Cassini 


range [m] 


31 


27324 


22 


Uranus 










Optical 


ra [arcsec] 


3849 


0.358 


0.351 


Optical 


de [arcsec] 


3835 


0.366 


0.361 


Neptune 










Optical 


ra [arcsec] 


3898 


0.368 


0.361 


Optical 


de [arcsec] 


3879 


0.360 


0.358 


Pluto 










Optical 


ra [arcsec] 


1023 


0.170 


0.170 


Optical 


de [arcsec] 


1023 


0.171 


0.171 



set. Observations deduced from the tracking of the Cassini 
spacecraft processed and provided by JPL (Folkner et al. 
2008) were also added to the set of observations. 

Global adjustments of planet initial conditions, Earth- 
Moon mass ratios, Astronomical Unit, Sun oblateness 
J2 as well as 34 asteroid masses were fitted to obtain 
INPOP08. Obtained values are presented in table [2] and 
table [3] as well as comparable values found in the liter- 
ature. The fit procedures are discussed in the following 
sections. 

5.1. Global results 

5.1.1. Contribution of VEX data 

The impact of VEX observations on the INPOP08 adjust- 
ment is very important. As one can see in Figure [3] and 
table [1] INPOP08 provides a much more accurate Venus 
orbit thanks to the VEX input data. This improvement is 
mainly induced by the fit of planet initial conditions and 
also by the improvement of the modclisation of the aster- 
oid perturbations on the inner planets. A mean dispersion 



of 4 meters (1-sigma formal dispersion) with a bias of 
1.5 meters is obtained after the fit, which corresponds 
to the level of accuracy of the VEX observations. This 
is an improvement in the estimation of the Earth- Venus 
distance of about a factor 42 compared to the previous 
ephemeris, INPOP06, which was not fitted to VEX data. 
The large effect that one can see on figure [3] is actually 
induced by the propagation of the previous planetary 
ephemerides uncertainties on Earth and Venus orbits. 
Such ephemerides as INPOP06, but also DE405, were 
only fitted to direct radar observations on the Venus 
surface (see Standish 1998b; Fienga et al. 2008) and to 
the VLBI data deduced from the tracking of the Venus 
mission Magellan in 1994. The Venus data sets were suf- 
fering from a big lack of observations since 1994 and from 
low accuracy of the direct radar observations. The com- 
bination of both explains the spectacular improvement 
observed on Figure [3l One can also notice in table [1] that 
even if the new Venus INPOP08 orbit induces significant 
changes in the values of VEX residuals compared to 
INPOP06 ones, the residuals obtained by comparison be- 
tween this new orbit and the observations used previously 
for the INPOP06 fit are quite similar to those obtained 
with INPOP06. This means that the VEX data are 
consistent with the old ones and their addition in the fit 
does not degrade the ephemerides on a large time interval. 



5.1.2. Mars observations and asteroids 

Using MEX tracking data and the INPOP06 Mars data 
set, we fitted the new asteroid ring modeling as well as the 
selection of asteroids described in section [31 Compared to 
INPOP06, the ring model was modified to a non-static 
ring, 3 more asteroids were included in the list of main 
perturbers and more asteroids have their masses fitted in- 
dividually in using a priori sigmas (Moyer 1971) of about 
30% constraints to their initial values. The choice of the 
fitted asteroid masses was done in a way to have only pos- 
itive masses even without the constraints. In INPOP08, 
34 asteroids have their masses fitted individually against 
5 in INPOP06. 

As one can see in table [T] and figure [H this new ap- 
proach and the input of information brought by the MEX 
data improve the ephemerides by reducing significantly 
the residuals. This improvement can be noticed with the 
new MEX data set which was not included in the pre- 
vious solution, but also with the MGS/MO data set in- 
cluded both in INPOP08 and INPOP06 adjustments. For 
the MGS/MO residuals we have an improvement of the 
results of a factor 4. 

Furthermore, as one can see on Figure [4] during the 
overlap period of the mission (from 2005.5 to 2006), the 
connection between MGS/MO and MEX data is done 
properly by INPOP08 with a bias of about 3.5 meters 
between MGS/MO and MEX. This bias is in the error es- 
timations presented in section [4] and the calibration done 
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by ESA before the MEX launch. Even if a great improve- 
ment can be noticed on Figure HI signals still remain in 
INPOP08 residuals due to a still possible lack in asteroid 
mass determinations but also to systematics induced by 
solar conjunction gaps in the data sets. 

In Table [2] and El asteroid masses fitted in INPOP08 
are compared to other published masses as well as val- 
ues for taxonomic densities and physical characteristics 
of the asteroid ring. In the case of INPOP08, the ring 
mass and distance were not fitted during the global ad- 
justment but during the construction of the ring (see sec- 
tion [3]) . Due to the direct correlation between the mass 
of the ring and its distance to the Sun, we decided to 
fix arbitrarily the distance at 3.14 AU. If this distance is 
chosen at 2.8 AU, the mass of the ring is estimated to 
be Mring — 0.3 ± 0.1 X 1O^^*'A/0 which is comparable 
with the GM estimated by DE414 and very close to the 
one given by INPOP06. In table [2 the taxonomic den- 
sities estimated by INPOP08 stay of the same order as 
those estimated by DE414 or DE421 and the values of 
masses estimated for the 5 biggest asteroids are also in 
good agreement with those published previously. Figure [5] 
depicts the same conclusion. Indeed, on figure [5l the val- 
ues estimated since 1991 (Baer and Chesley 2008) of GMs 
of Ceres, Pallas and Vesta are plotted as well as the val- 
ues provided by INP0PG6, DE414, DE421 and INPOP08 
which correspond to the four last points dated in 2008 
and 2009. On this graph, are represented (using darked 
circles) asteroid masses estimated by close-encounters (see 
for instance Baer and Chesley 2008, Hilton 1999), asteroid 



masses obtained after radar imaging (Michalak 2000) and 
masses fitted during the global adjustments of planetary 
ephemerides marked as black crosses. The close-encounter 
method was used a lot in the 90s but due to the difficulty 
of the method (accurate observations of the perturbed ob- 
jects to obtain very accurate estimation of its orbit and 
to survey the orbit during the close-encounters), the ef- 
fort was not maintained during the last decade. On the 
other hand, the evolution of planetary ephemerides gives 
very stable estimations for the largest asteroids with un- 
certainties smaller than the close-encounters determina- 
tions. In 2008, (Baer and Chesley) published values es- 
tablished on the basis of several encounters and obtained 
values of masses for Ceres, Pallas and Vesta very close to 
the ones obtained by the planetary ephemerides. For Ceres 
and Vesta, the uncertainties of the close-ecounters deter- 
minations are at the same level of accuracy as planetary 
ephemerides ones. INPOP08 gives original values offset 
from the previous estimations but at the limit of the error 
bars. 

In table [3l are given the masses of the 28 asteroid, 
other than the 5 bigs given in table [H fitted individually 
in INPOP08. The estimated mass values are put in the 
table as their mean impact on the Earth-Mars distances 
(estimated over 20 years) decreases. For some of them, the 
INPOP08 masses agree in the limit of 30 % with values 
estimated previously either by close-encounters methods 
(as for (16) Psyche, (52) Europa, (88) Thisbe, (8) Flora, 
(15) Eunomia, (18) Melpomene (19) Fortuna, and (21) 
Lutecia) or by adjustment of planetary ephemerides (as 
for (41) Daphne, (29) Amphitrite, (409) Aspasia, (704) 
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Fig. 4. MEX and MGS/MO 1-way residuals with the new INPOP08 fitted to VEX, MEX and MGS/MO data (right- 
hand side plots) and INPOP06 fitted to all MGS/MO and a part of MEX data (left-hand side plots). The improvement 
of INPOP08 is mainly induced by changes in the asteroid mass determinations and the addition of MEX data, see 
section |3] in the text. The y-axis unit is meter. 



Interamnia, or (532) Herculina). However, for some oth- 
ers, the differences from previous estimations can reach 
500 % as for the case of (9) Metis or even worse (as for (6) 
Hebe) . We estimate to 50% of estimations close to previ- 
ously published masses (differences below 50 %), 21% of 
first estimations and 4 very unsufficent estimations (dif- 
ferences greater than 500 %). 

Such differences can be explained by the differences 
in the modelisation used to fit the data. For instance, by 
the differences in the data reduction, the weight used in 
the fit as well as by fitting individually a different number 
of asteroids and by absorbing the effect of all the other 
minor bodies by the fit of the taxonomic classes and the 
addition of a ring, one absorbs in the fit of individual 
masses dynamical effects which are not really caused by 
one particular asteroid but by a group of asteroids induc- 
ing effects very similar in amplitudes and periods. The 
fitted mass for this particular asteroid represents then not 
only its gravitational potential but also the one induced 
by other asteroids. This remark can be illustrated by the 
unrealistic values of some densities deduced from masses 
and presented in table [3] as well as the almost zero mass of 
the asteroid (747) Winchester, the perturbations induced 
by (747) being quite negligible or not clearly separated 
from another asteroid perturbations like in the case of (6) 
Hebe. In order to better estimate the real mass value of 
one particular asteroid, one could isolate the arc of Mars 
orbit where and when this asteroid has the biggest impact 



and try to estimate its mass only on that arc. Such an 
investigation is presented in (Somenzi et al. 2009). 

5.1.3. Cassini normal points 

Cassini normal points for Saturn were provided by JPL 
and are described in details by (Folkner et al. 2008). The 
impact of these data is not negligeable. They give very im- 
portant informations about the Saturn orbit itself by pro- 
viding very accurate estimations of the Earth-Saturn dis- 
tances and geocentric angular positions of Saturn. Thanks 
to these, our knowledge of this orbit especially in term of 
distances from the Earth reduces from about 27 km to 22 
meters at the epoch of Cassini normal points. 

5.2. AU fixed and reseated GM of the sun version 

As a first step to a new generation of planetary 
ephemerides fitting the GM of the sun and using a fixed 
value of AU, an alternative version of INPOP08 was built. 
This version is based of the INPOP08 ephemeris presented 
in the previous section. However, instead of providing a 
value of AU fitted to observations, the AU is fixed to the 
(lERS 2003) value and the change is absorbed by the mul- 
tiplication of all the initial conditions of planets and as- 
teroids as well as all the masses including the GM of the 
sun by a factor equivalent to AUfittod over AUiersos- A 
new value for the GM of the Sun is then deduced and pre- 
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Table 2. Physical parameters fitted in INPOP08. Other values deduced from planetary ephemerides are presented for 
comparison. EPM08 stands for (Pitjeva 2008). The uncertainties are given at 5-sigma for the GMs and 1 formal sigma 
for others, n/a stands for non-avalaible and NE as non-estimated. For AU, the presented values are the differences in 
meters between the fitted values and the AU value of the lERS conventions 2003, AUiersos = 149597870.691 km. 



Unit DE414 DE421 EPM 2008 INPOP06 INPOP08 



Mass of Ceres 


10" 


10 


o 


4 


699 ± 0.028 


4.685 


4.712 ± 0.006 


4.756 ± 0.020 


4 


658 ± 0.045 


Mass of Vesta 


iU 


10 




1 


358 ± 0.016 


1.328 


1.344 ± 0.003 


i.o48 ± U.UiO 


1 


392 ± 0.015 


Mass of Pallas 


10" 


10 


Mq 


1 


026 ± 0.028 


1.010 


1.027 ± 0.007 


1.025 ± 0.005 


1 


076 ± 0.010 


Mass of Juno 


10" 


10 


Mq 





149 ± 0.015 


0.116 


n/a 


NE 





075 ± 0.015 


Mass of Iris 


10" 


10 







060 ± 0.010 


0.060 


n/a 


0.058 ± 0.005 





050 ± 0.010 


Mass of Bamberga 


10" 


10 







047 ± 0.007 


0.048 


n/a 


0.046 ± 0.003 





056 ± 0.004 


Mass of Ring 


10" 


10 







30 


NE 


n/a 


0.34 ± 0.1 


1 


± 0.3 


Distance of Ring 


UA 






2 


8 


NE 


n/a 


2.8 


3 


14 


Density of the C class 








1 


62 ± 0.07 


1.09 


n/a 


1.56 ± 0.02 


1 


54 ± 0.07 


Density of the S class 








2 


08 ± 0.19 


3.45 


n/a 


2.18 ± 0.04 


1 


94 ± 0.14 


Density of the M class 








4 


32 ± 0.37 


4.22 


n/a 


4.26 ± 0.12 


4 


98 ± 0.50 


Sun J2 


lo- 


7 




2 


34 ± 0.49 


2.0 


n/a 


2.46 ± 0.40 


1 


82 ± 0.47 


EMRAT 








81.300568 


81.300569 81.3005690 ± 0.0000001 


NE 


81.300540 ± 


AU-AUiersos 


rn 






9 


8 ± 0.15 


8.6 ± 0.15 


4.4 ± 0.10 


NE 


8 


22 ± 0.11 



_ 1 


1 ' 1 ' 1 ' 1 ' 


' 1 ' 1 

Ceres 


" , 1 , 


'I f' — ^ 

1 . 1 < 1 . 1 . 


,1,1," 




Table 4. AU and GM of the sun values used for the 
construction of INPOP08 and INPOP08b. In the case of 
INPOP08, the AU is fitted and the mass of the Sun fixed 
to the DE405 one. In the case of INPOP08b, the AU is 
fixed to be equal to AUiersoS: a new estimation of 
the mass of the Sun is deduced as described in section [5T2l 




1998 2000 2002 2004 2006 2008 2011 

years of publication 

Fig. 5. Estimations of Ceres, Pallas and Vesta masses 
found in the literature. The GMs (in unit of lO"^" so- 
lar mass y-axis) are plotted by years of publications (x- 
axis) . The dark red circles represent the GMs estimated in 
using close-encounters as described by Baer and Chesley 
(2008). The crosses represent GMs fitted with the plan- 
etary ephemerides. The upper triangles represent masses 
deduced from radar observations as described by Michalak 
(2000). The last crosses corresponding to 2009 give values 
obtained by INPOP08. 



sented on Table H This version of INPOP08 is equivalent 
to the previous one but allows the next versions of plane- 
tary ephemerides to fit directly the GM of the sun instead 
of the AU. 



AU 
km 



GM0 
km'\s"2 



Fitted DE405 
149597870.69922 132712440039.87900 



INPOP08 



IERS03 Deduced 
INPOP08b 149597870.69100 132712440039.878750 



6. Estimation of uncertainties 

Often, users are asking for the accuracy of planet 
ephemerides. There are two methods to answer. The first 
one is to compare the ephemerides to observed positions 
not used in the fit, and preferably to observations situated 
out of the time interval of observations used for the adjust- 
ment. Such results give information to the users who di- 
rectly deal with observations of planets or satellites: what 
accuracy in millarcseconds can we predict the position of 
such a planet in right ascension and declination, what ac- 
curacy in meters obtained for the estimation of Earth- 
Mars distances etc... . In section l6Tl we present such an 
extrapolation of INPOP08 compared to observations not 
used in the fit and usually outside the fit interval. 
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Table 3. Asteroid GMs fitted in INPOP08. Other values deduced from planetary ephemerides are presented for 
comparison. Reference [1] stands for (Baer and Chelsey, 2008). [2] are values estimated in DE414 but not published 
with crrorbars. They arc directly extracted from the DE414 header. [3] stands for (Bacr ct al. 2008) and [4] for DE421 
values as published by (Folkner et al. 2008). The given uncertainties are formal accuracy given at 1-sigma. Column 4 
and 5 are dedicated to physical characteristics of the objects: in column 4 are given radius in kilometers, and spectral 
type in column 5. In column 6, arc given the densities obtained in using the given radius and the estimated masses 
given in column 2. Density errorbars are estimated with the given mass 1-sigma uncertainties and accuracies on radius 
when provided by the PDS website. The last column gives the maximum impact of each asteroid on the Earth-Mars 
distances over the 1998 to 2008 intervalle of time. 





TNPDPns 




r 


type 


p 






1 n~i^ 








_3 

g.cm 


[ml 
[mj 


16 Psyche 


1.596 ± 0.032 


1.29 ± 0.17 [1] 


126.6 ± 2.0 


M 


3.7 ± 0.1 


138 






1.688 [4] 










29 Amphitrite 


0.491 ± 0.09 


1.00 ± 0.35 [1] 


106.1 ± 3.4 


S 


1.9 ± 0.4 


118 






0.684 [4] 










14 Irene 


0.071 ± 0.01 


0.2623 [2] 


76.0 ± 8.0 


S 


0.8 ± 0.3 


69 






0.413 ± 0.073 [3] 














0.263 [4] 










704 Interamnia 


1.623 ± 0.01 


3.58 ± 0.42 [1] 


158.3 ± 2.0 


c 


1.9 ± 0.1 


65 






1.862 [4] 










532 Herculina 


0.546 ± 0.01 


0.667 [2] 


111.2 ± 2.0 


s 


1.9 ± 0.1 


63 






0.669 [4] 










6 Hebe 


0.016 ± 0.011 


0.759 ± 0.142 [1] 


92.5 ± 1.5 


s 


0.1 ± 0.1 


63 






0.255 [2] 














0.457 [4] 










52 Europa 


1.724 ± 0.100 


0.976 ± 0.22 [1] 


151.25 ± 2.7 


c 


2.4 ± 0.2 


52 






1.023 [4] 










9 Metis 


0.115 ± 0.100 


1.03 ± 0.24 [1] 


95.0 ± 9.5 


s 


0.6 ± 0.5 


49 






0.600 [2] 














0.428 [4] 










19 Fortuna 


0.202 ± 0.02 


0.541 ± 0.008 [1] 


101.7 


c 


0.9 ± 0.1 


48 






0.350 [4] 










128 Nemesis 


0.169 ± 0.033 




94.1 ± 2.0 


c 


1.0 ± 0.2 


42 


18 Melpomene 


0.091 ± 0.02 


0.151 ± 0.051 [3] 


70.3 ± 1.5 


s 


1.2 ± 0.3 


40 






0.201 [4] 










15 Eunomia 


2.230 ± 0.018 


1.68 ± 0.08 [1] 


127.7 ± 7.0 


s 


5.0 ± 0.7 


37 






1.239 [4] 










88 Thisbe 


0.863 ± 0.026 


0.57 ± 0.18 [1] 


116.0 ± 


c 


2.6 ± 0.1 


36 


409 Aspasia 


0.105 ± 0.003 


0.163 [4] 


80.8 ± 4.0 


c 


0.9 ± 0.1 


29 


216 Kleopatra 


0.353 ± 0.012 


0.255 [2] 


67.5 ± 


M 


5.3 ± 0.2 


27 






0.226 [4] 










21 Lutetia 


0.1034 ± 0.03 


0.129 ± 0.012 [3] 


47.9 ± 4.1 


M 


4.4 ± 1.3 


24 






0.105 [4] 










31 Euphrosyne 


2.99 ± 0.68 


1.54 [2] 


127.95 ± 5.5 


c 


6.7 ± 1.6 


23 






0.313 ± 0.059 [3] 














0.860 [4] 










23 Thalia 


0.003 ± 0.001 




53.8 ± 2.2 


s 


1.06 ± 0.08 


23 






0.097 [4] 










354 Eleonora 


0.488 ± 0.035 


0.246 [2] 


77.6 ± 4.2 


s 


4.9 ± 0.4 


22 






0.247 [4] 










192 Nausikaa 


0.137 ± 0.041 


0.081 [4] 


51.6 ± 1.0 


s 


4.6 ± 1.4 


22 


747 Winchester 


0.0004 ± 0.0002 


0.148 [4] 


85.9 ± 1.5 


c 


3e-3 ± lc-3 


22 


129 Antigone 


0.717 ± 0.014 




104.58 ± 


M 


3.0 ± 0.1 


21 


41 Daphne 


0.527 ± 0.05 


0.467 [2] 


87.00 ± 6.0 


c 


3.8 ± 0.7 


20 






0.398 [4] 










173 Ino 


0.366 ± 0.030 




77.0 ± 1.6 


c 


3.8 ± 0.3 


20 


89 Julia 


0.359 ± 0.014 




75.8 ± 1.5 


c 


4.0 ± 0.4 


15 


8 Flora 


0.535 ± 0.005 


0.178 [2] 


67.9 ± 1.0 


s 


8.1 ± 0.5 


14 






0.426 ± 0.045 [3] 














0.178 [4] 










12 Victoria 


0.117 ± 0.003 




56.4 ± 1.5 


s 


3.1 ± 0.2 


13 


139 Juewa 


0.359 ± 0.010 


0.142[4] 


78.3 ± 1.5 


c 


3.6 ± 0.10 


10 
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Fig. 6. Residuals obtained by comparisons between INPOP08 and MEX tracking data (left hand side) and VEX 
(right hand side) observations. The small dark points represent residuals of observations which were used in the fit 
of INPOP08 and the big dark points represent residuals of observations not used in the fit. The large slope of MEX 
residuals may possibily be improved in the future with better determination of asteroid masses. 



However, for all the other users, the comparison to ob- 
servations is interesting only if they correspond to their 
specific use. In general, users prefer to have a more 
theoretical estimation such as the uncertainties in the 
barycentric positions and velocities of the Earth. To an- 
swer this question, we use comparisons between different 
ephemerides to estimate the level of variations induced by 
an improvement of modeling or adjustment (section 16. 2|) . 



6.1. Extrapolation tests 

Two sets of observations of MEX and VEX were kept out 
from the fit procedure as test sample of INPOP08. These 
sets were obtained from April 2008 to September 2008. 
Differences in meters obtained by comparisons between 
the observed and INPOP08 distances can be found in fig- 
ure [6l For Mars, the linear drift observed is in the order 
of error expected (about 20 meters over half a year) and 
is mainly due to asteroid modeling and adjustment. 

For Venus, the distribution is noisier, but in good 
continuity with the postfit residuals obtained after the 
INPOP08 fit and presented in Figure [6l The biggest no- 
ticeable effect on these residuals is induced by the VEX 
solar conjunction that occured on June 9th 2008. This so- 
lar conjunction has introduced unmodeled noise into the 
doppler and ranging data of VEX during a period of few 
days (between the 06 to 11 June inclusive) during when 
the spacecraft tracking was stopped. The peak that can 
be noticed on Figure [6] is then mainly induced by this 
conjunction. 



6.2. Uncertainties by comparisons with other 
ephemerides 

Another way to estimate the uncertainties of planet or- 
bit is to compare several ephemerides. By such compar- 
isons, one estimates first the impact of the differences 
in dynamical modeling of each ephemeris and in adjust- 
ment to observations: the ephemerides could have been 
fitted over different sets of observations with different sets 
of weight and with different sets of parameters. Table [5] 
presents the maximum differences in barycentric positions 
and velocities between INPOP08 and INPOP06. Tabled 
presents the corresponding maximum differences between 
INPOP08 and DE421. Two intervals of time are consid- 
ered to estimate these differences: interval A from 1990 
to 2010 and interval B from 1900 to 2050. Furthermore, 
plots of heliocentric differences in ecliptic longitudes ob- 
tained by comparisons between the same ephemerides are 
presented in figure [71 The use of Cassini normal points in 
DE421 and INPOP08 make the differences between the 
two ephemerides decrease significantly for the longitude 
of Saturn, as well as for the other outer planets. 

If one consider that these differences give a good es- 
timation of their actual accuracies, then thanks to the 
Cassini data, we can see the improvement of the orbit at 
the epoch of these observations. 

Limited improvements of the outer planet orbits will 
be possible by the addition of more accurate data in the 
fit, mainly tracking data of spacecraft orbiting or cruising 
one of these systems, such as the New Millenium mission. 
Important improvements will be obtained after having sig- 
nificant portion of orbit covered by accurate observations 



Fienga et al: INPOP08, a 4-D planetary ephemeris 



13 



Table 5. Maximum differences between INPOP06 and 
INPOP08 estimated on interval A from 1990 to 2010 and 
on interval B from 1900 to 2050. 



Barycentric Barycentric 
Positions (m) Velocities (m/d) 



Interval 


A 


B 


A 


B 


Mercury 


3500 


4600 


260 


280 


Venus 


1700 


2500 


37 


41 


EMB 


760 


4000 


11 


48 


Mars 


980 


14000 


9 


140 


Jupiter 


240000 


1300000 


350 


1900 


Saturn 


320000 


560000 


161 


320 


Uranus 


1800000 


1800000 


290 


410 


Neptune 


2400000 


7500000 


320 


620 


Pluto 


7500000 


140000000 


1700 


7200 



Table 6. Maximum differences between DE421 and 
INPOP08 estimated on interval A from 1990 to 2010 and 
on interval B from 1900 to 2050. 



Barycentric Barycentric 
Positions (m) Velocities (m/d) 



Interval 


A 


B 


A 


B 


Mercury 


3900 


4400 


290 


310 


Venus 


740 


3200 


13 


39 


EMB 


870 


3900 


11 


36 


Mars 


1300 


9200 


10 


82 


Jupiter 


210000 


1300000 


290 


2000 


Saturn 


35000 


40000 


22 


23 


Uranus 


1400000 


2500000 


190 


390 


Neptune 


3700000 


14000000 


480 


1200 


Pluto 


4300000 


120000000 


830 


6700 



which will require at least five years of continuous obser- 
vations for Jupiter orbit. 

For Mars and the EMB, one can see on figure[7]the dif- 
ferences are also reduced between DE421 and INPOP08 
compared to the differences obtained between INPOP06 
and INPOP08. As well for the EMB in table 5, if one uses 
the differences between ephemerides as a tool to estimate 
ephemerides accuracy, the uncertainties on its barycentric 
positions are about 900 meters over 10 years and 4 kilo- 
meters over 100 years, and about 10 meters per day to 50 
meters per day for its barycentric velocities. 

7. Perspectives and conclusions 

In this article, INPOP08 was fully described as a new 4-D 
planetary ephemeris. 

New modeling of asteroid perturbations were intro- 
duced and as described by Kuchynka et al. (2009), im- 
provements in the methods of asteroid mass determina- 
tions are possible. Introduction of these new methods will 
be done in the next INPOP version, leading to an improve- 
ment of the extrapolation capabilities of the ephemerides 
as well as more realistic determination of asteroid masses. 



New sets of observations were used for the construction 
of INPOP08: the Mars Express and Venus Express track- 
ing data processed by ESA and normal points deduced 
from the Cassini mission. Thanks to these data, large im- 
provements were made in the determination of the orbits 
of Venus and Saturn. They open doors to new tests of 
gravity which were hitherto limited to only Earth-Mars 
distance accurate observations. These tests and their re- 
sults will be presented in a dedicated paper very soon. 

However, as was already stressed in section [5. 1. 11 VEX 
tracking data and Saturn Cassini normal points have a 
crucial role to play in such investigations. 

The LLR data have been used for the adjustment of 
the lunar orbit and libration in INPOP08 (Manche et al. 
2007). The fit to these data as well as the complete lunar 
model will be described in a forthcoming paper. 

The INPOP08 planetary ephemeris is available for 
users on the INPOP website (www.imcce.fr/inpop) as 
Chebychev polynomials coefficients interpolated via pro- 
vided routines in C and fortran. Positions and velocities 
of the nine planets, the Sun, the moon libration and the 
3 Euler angles of the Earth rotation can be extracted as 
well as the TT-TDB relation at any time. An alternative 
version with rescaled GM of the sun and AU fixed is also 
provided on request. 

In the next version, we plan to adjust directly the GM 
of the sun instead of the AU. Simulations of Mercury mis- 
sions data and use of fly- by data will be introduce in order 
to evaluate the improvements in terms of accuracy limits 
in the determination of PPN parameters and J2. 
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